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Abstract — In this paper, we propose three novel modulation 
techniques, i.e., concentration-based, molecular-type-based, and 
molecular-ratio-based, using isomers as messenger molecules 
for nano communication networks via diffusion. To evaluate 
achievable rate performance, we compare the proposed tech- 
niques with conventional insulin based concepts under practical 
scenarios. Analytical and numerical results confirm that the 
proposed modulation techniques using isomers achieve higher 
data transmission rate performance (max 7.5 dB signal-to-noise 
ratio gain) than the insulin based concepts. We also investigate 
the tradeoff between messenger sizes and modulation orders 
and provide guidelines for selecting from among several possible 
candidates. 

Index Terms — Nano communication networks, molecular com- 
munication, modulation technique, isomer, diffusion, and mes- 
senger molecule. 



I. Introduction 

SINCE Richard Feynman's talk on top-down nanotech- 
nology, people have eagerly pursued practical work on 
smaller and smaller scales, most notably through nanotech- 
nology |2|. Nanotechnology, recently, produced a new branch 
of research called nano communication networks (NCNs) |[3J- 
(51. NCNs interconnect several nanoscale machines (a.k.a. 
nanomachines) to perform simple tasks or to carry out more 
complex tasks in a cooperative manner |3|, (4j. These net- 
works are not just smaller versions of traditional communica- 
tion networks; they have their own features and are applicable 
in many fields, including biomedical, industrial, military, and 
environmental fSl. 

NCNs can be realized by several methods. We can rely on 
traditional communication systems that use electromagnetic 
fields or ultrasonic waves. Such a method, however, has to 
overcome some radio frequency (RF) device barriers |[3|. On 
the practical front, researchers are considering new materials, 
such as carbon-nano tubes (CNTs) and graphene |[6|, ||7|. 
Areas that have yet to receive much research, however, are 
channel models and human body absorption with terahertz 
bands. Therefore, researchers have introduced a new concept 
utilizing diffusion that is especially useful for short range 
communication |[8J. 
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This new paradigm of communication, molecular communi- 
cation, sends/receives information-encoded molecules between 
nanomachines |[3|, |[9l-|T2|. Several papers have addressed the 
achievable transmission rates (achievable rate hereafter) of the 
communication system, and the author in |4| evaluated, using 
a circuit model, the normalized gain and delay of the system. 
Though a great deal of these papers have ignored possible 
noise sources, recently several works have focused on noise 
analysis |[13)-|[T5). 

One of the advantages of molecular communication is its 
biocompatibility. By using bio-molecules, we can create inher- 
ently biocompatible systems that require no harmful inorganic 
materials. Indeed, biocompatibility is regarded by many as the 
most important and challenging issue concerning intra-body 
applications. Molecular communication, driven by chemical 
reactions, is also energy efficient 1 12|. For these reasons, this 
paper mainly focuses on molecular communication. 

The authors in p6| , |T7| studied extensively the funda- 
mentals of molecular communication via diffusion. In |16|, 
they investigated a new energy model to understand how 
much energy is required to transmit messenger molecules 
and 1 17] introduced concentration-based and molecular-type- 
based modulation techniques. The researchers also have com- 
pared, by using a simple binary symmetric channel model, 
the achievable rate. In analyzing their modulation techniques, 
however, |T6| , p7| did not clearly suggest concrete structures 
for the messenger molecules. They considered insulin-based 
nano networks, but it is unclear how these could be utilized in 
practice. In this paper, so as to maximize the achievable rate 
with less transmit power/energy, we propose using isomers 
as messenger molecules. We also suggest a new ratio-based 
modulation method. To consider the properties of isomers, we 
slightly modify the energy model described in 1 16 |. To the best 
of our knowledge, the proposed method is the first attempt to 
design appropriate practical messenger molecules for NCNs 
via diffusion. 

This paper is organized as follows. Section [ll| describes 
the channel and energy model under consideration. Sec- 
tion Inl] explains isomers for messenger molecules and also 



proposes three modulation techniques, i) concentration-based, 
ii) molecular-type-based, and iii) ratio-based. We present the 
numerical results in Section |IV| and Section |V| offers our 
conclusions. 
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Fig. 1. System model of nano communication via diffusion with a single 
transmitter and a single receiver. 



II. System Model 

We consider a nano communication system that, for sim- 
plicity, consists of a single transmitter and a single receiver, 
as illustrated in Fig.[T] The information-encoded molecules are 
called messenger molecules, and they can propagate through 
the medium by several methods. Three main NCN archi- 
tectures are based on the propagation type (e.g., walkaway- 
based, flow-based, and diffusion-based) |4|. For this system, 
the diffusion-based propagation is chosen for analysis, so the 
messenger molecules diffuse through the medium (i.e., liquid) 
at body temperature. This paper assumes that no collisions 
occur among the propagating messenger molecules. 

A. Channel Model 

The particles or messenger molecules released from a 
transmitter nanomachine spread out through the medium by 
Brownian motion p3| , fT6| . Such motion is basically driven 
by diffusion, meaning the particles move from areas of higher 
concentration to areas of lower concentration, and the displace- 
ment of messenger molecules follows a normal distribution 
with zero mean. Fig. [2] is one example of one-dimensional 
Brownian motion. Also, Fig. [3] shows molecular signal attenu- 
ation as functions of distance and time. The standard deviation, 
cr, of the displacement can be obtained as follows: 




-x'^ /ADt 



(1) 
(2) 

(3) 



where c indicates the concentration of Brownian particles at 
time t at point x, and D represents the diffusion coefficient 
of the particles calculated from the following: Boltzmann con- 
stant (Kb), the temperature (T), the viscosity of the medium 
(77), and the radius of a messenger molecule (r^^). Eq. ([T]) 
is Fick's second law of diffusion. By solving this partial 
differential equation, we obtain the general equation for the 
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Fig. 2. An example of the Brownian motion of a particle in \-D motion. 

concentration of particles represented in ([2]). The first moment 
of ([2]) is zero, which indicates the displacement has a zero 
mean, and the second moment becomes the variance. Hence, 
the displacement has a standard deviation of V2Di, and finally 
we have 



AX - A/'(0,2Dt). 



(4) 



When n number of messenger molecules are transmitted by 
the transmitter nanomachine, the molecules have a probability 
of hitting the receiver nanomachine. We represent this as 
a binomial distribution of n times of trials with a hitting 
probability (Phit) for each trial. Phit is determined by the 
symbol duration (Tg) and the distance between the transmitter 
and the receiver (d), which are both affected by the diffusion 
coefficient (D). If n is large enough, and nPhu is not zero, 
the binomial can be approximated as a normal distribution. In 
addition, we have to take into account the overflow molecules 
from the previous symbol as well as the molecules sent from 
the current symbol]^ Thus, the total number of molecules 
during one symbol duration differs according to the current 
and previous symbols. The following represents the number 
of molecules from the current and previous symbols and the 
noise term. 

Binomial{n^ Phit{d^ Ts)) 

- M{nPM Ts). nPu,{d, Ts){l - PM Ts))). 
Nc^JV{npunpi{l-pi)), (5) 
Nr, J\f{np2,np2{l -P2)) - M{npi,npi{l -pi)), 
^A/'(0,a2), 



Nn 



where, Nc denotes the number of molecules transmitted and 
received during the current symbol duration, and Np denotes 
the number of molecules transmitted in the previous symbol 
duration but received during the current symbol duration. For 

^In our paper, we consider the overflow only from the previous symbol since 
Phit for 2Ts shows a significant increase compared to Phit for Ts , but, as more 
time passes, no meaningful difference exists(i.e., memory channel) 1 16|. 
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Fig. 3. Signal attenuation as functions of distance and time. 

simplicity, we use pi for Phit during and p2 for Phit during 
2T5. The molecules from noise sources can be summed up 
as Nn, and, in this paper, it is assumed as Additive White 
Gaussian Noise (AWGN) |[17}. 

B. Energy Model 

We use the energy model described in fT6l and assume 
that basic eukaryotic cells are a perfect model for nanoma- 
chines. First, the messenger molecules are synthesized inside 
a nucleus and encapsulated in vesicles during intra-cellular 
propagation. The vesicles are then carried to the boundary of 
the transmitter nanomachine. Lastly, the messenger molecules 
are released into the propagation medium through membrane 
fusion of vesicles. Each step is illustrated in Fig. |4] Fig. |4] 
represents a simplified form of an eukaryotic cell structure 
to model a nanomachine. The large circle represents the cell 
membrane filled with the medium (e.g., cytoplasm), and the 
inner circle represents the nucleus of a cell. The structure 
around the nucleus is a Golgi apparatus that produces and 
transports vesicles. Furthermore, some interconnections exist 
between the Golgi and the cell boundary-microtubules or 
molecular rails-which messenger molecule-carrying vesicles 
travel along. In the cell membrane, several gated channels 
exist that the nanomachine can, as needed, control the on/off 
of them. Specifically, the gated ionic channel can be on/off 
depending on the presence of ions near the channel. 

The energy needed in each step is calculated. Moreover, 
we can obtain the cost of synthesizing messenger molecules 
(hexoses as an example) by the enthalpy of formation (AH, 
which is the change of enthalpy that accompanies the forma- 
tion of one mole of the compound from its elements). The 
final energy model can then be presented as follows: (see |T6| 
for more details) 

TL 

Et = nEs + —{Ev ^Ec^ Ee), 

Cy 

7-, AH}iexose T 1 1 

Es = 7^ J per messenger molecule, 




Fig. 4. Energy model of nano communication composed of four steps. The 
first step (1) is for generating the messenger molecules, (2) for synthesizing 
the vesicles, (3) is for carrying them, and (4) is for extracting them to the 
medium. 
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where, Et is the total energy cost required to transmit n 
number of molecules, Es is the synthesizing cost of one 
hexose molecule calculated from the sum of bond energies 
(e.g., the enthalpy change), and Ey is the vesicle-synthesizing 
cost having a radius of Ty. Ec is the cost of intra-cellular 
transportation having cell radius of runit^ Ee is for membrane 
fusion, and Cy is the vesicle capacity, the number of messenger 
molecules one vesicle can carry, which is related to the radius 
of messenger molecules rmm- 

III. Modulation Techniques 

In molecular communication via diffusion, the unique prop- 
erties of the messenger molecules can determine modulation 
techniques. Two such techniques proposed in |17| include 
the use of concentration and types of messenger molecules. 
A usable messenger molecule itself, however, has yet to 
be proposed. Thus, this paper proposes practical messenger 
molecules as well as a new modulation technique using ratios 
of the molecules. We also analyze and compare achievable 
rates by applying different modulation techniques. 

A. Isomers for Messenger Molecules 

The most important thing when designing messenger 
molecules is that they have to be non-toxic to the human body. 
The one suggested in 1 16], however, is highly flammable (i.e., 
hydrofluorocarbon), which means it may be inappropriate for 
practical applications. 

For several reasons, potential candidates for messenger 
molecules are isomers, molecules composed of the same 
number and types of atoms 1 18 1. First of all, they consist of the 
same type of atoms, lightening the burden at the transmitter 
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H^^^O H^^^O H^^^O H^^^O H^^^O H^^^O H^^^O H^^^O 

I I I I I I I I 

H-C-OH HO-C-H H-C-OH HO-C-H H-C-OH HO-C-H H-C-OH H-C-OH 

I I I I I I I I 

H-C-OH H-C-OH HO-C-H HO-C-H H-C-OH H-C-OH HO-C-H HO-C-H 

I I I I I I I I 

HO-C-H HO-C-H H-C-OH H-C-OH H-C-OH H-C-OH HO-C-H HO-C-H 

I I I I I I I I 

H-C-OH H-C-OH H-C-OH H-C-OH H-C-OH H-C-OH H-C-OH H-C-OH 

I I I I I I I I 

CH2OH CH2OH CH2OH CH2OH CH2OH CH2OH CH2OH CH2OH 

D-Allose D-Altrose D-Glucose D-Mannose D-Gulose D-ldose D-Galactose D-Talose 

Fig. 5. D-form isomers of aldohexoses. 
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Fig. 6. Mirror image of D- and L- glucose. They are enantiomers. 
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nanomachine that synthesizes them. For numerical analysis, 
this paper uses the isomers known as hexoses, especially 
aldohexoses (i.e., aldose forms of the hexoses). Note that we 
can select one from the aldoses family (e.g., hexoses, pentoses, 
tetroses, or trioses) based on the required modulation order. 

Aldohexoses are monosaccharides with the chemical for- 
mula of CqHi20q. They have four chiral carbon atoms, 
not superposable on the mirror images, which give them 16 
(=2^) stereoisomers [19| . Fig. [s] illustrates eight kinds of 
D- form diastereomers. The enantiomers of each molecule 
are another set of eight L- form diastereomers. Therefore, 
aldohexoses have 16 different shapes. Here, diastereomers are 
stereoisomers that are not enantiomers (mirror-image isomers). 
For example, I^-glucose and L-glucose as shown in Fig. |6] are 
enantiomers. On the other hand, I^-glucose and I)-galactose- 
isomers but not mirror images-are diastereomers. 

When each isomer is dissolved in an aqueous solution, it 
mostly exists as a cyclic form (i.e., a ring form). I^-glucose, for 
instance, undergoes nucleophilic addition reaction generating 
four cyclic anomers, a-, j3- forms of D- glucopyranose and D- 
glucofuranose |20|. We consider, however, only two pyranose 
(a-, P-) forms since they predominate in the solution with 
36 and 64 percentages, respectively. If the functional groups 
attached to carbon number 1 (CI) and C5 shown in Fig.[7]have 
a trans-structure, it is called a- form, and if a cis-structure, 
form. They interconvert each other in solution through a 
process called mutarotation. Therefore, by deploying hexoses 
groups into the system, there are 32 different isomers in total. 



Fig. 7. D -glucose and its a- and /3- anomers. They undergo rapid 
interconversion in water, which is called mutarotation. 



B. Molecular-Concentration-Based 

When the concentration of messenger molecules is used to 
differentiate symbols, the technique is known as concentration- 
based modulation, originally introduced in |17|. Different 
concentration values represent different symbols, which are 
divided by proper thresholds. This is called concentration 
shift keying (CSK), and also can be referred to as isomer- 
based CSK (ICSK) here since we use isomers as the mes- 
senger molecules. In binary-ICSK (B-ICSK), one threshold 
exists to represent two different symbols, and a receiver 
nanomachine decodes the symbol as '1' if the number of 
the received messenger molecules exceeds the threshold (r), 
'0' otherwise. Generally, 2"^ -ICSK systems transmit n bits 
per symbol, and requires 2^-1 number of thresholds. In this 
system, there is, theoretically, no limit in the modulation 
order. As the modulation order increases, however, so does 
the error probability since the minimum distance between two 
neighboring thresholds decreases. The system only uses one 
kind of molecule, and I^-glucopyranose is used for analysis 
here. 

1) B-ICSK: 



Pa(0,0)=P5(0,0,0)+P5(l,0,0) 
= (^) [P{Nn<T)^P{N^^Nn<T) 



KIM AND CHAE: NANO COMMUNICATION NETWORKS VIA DIFFUSION 



5 



r - n{p2 -pi) 



+ 1-Q, , 

Pa(0,l) =Pb(0,0,l) + Pb(l,0,l) 
= (2) [PiNn >t) + P{Np + N^>T 



(6) 



Pa(l,0) = Pfc(0,l,0) + P6(l,l,0) 



= (^) [p{N, + N„<t) + P{Nj 



Ne + N„< t) 



T — npi 



1-Q 



V"[P2(1 - P2) + 2pi (1 - Pi)] + ct2 ^ 

Pa(l,l) =P6(1,1,1) + P6(0,1,1) 



■ iVe + iV„ > r) + P(Are + A/"™ > r) 



Q 



np2 



^ Vn[p2 (1 - P2) + 2pi (1 - pi)] + a2 ^ 
- - npi 



where, Pa(-^, ^) indicates the probabiHty of X sent and Y 
received, and P5(Z, X, F) is the probability of X sent, Y 
received, and Z previously sent. Q{-) is the tail probability 
of the normal distribution. In ([6]), (a) results from ^ and the 
relationship between Q function and the normal distribution. 

2) Q-ICSK: In Q-ICSK systems, three threshold values (ri, 
r2, and r^) are required to represent four different symbols. 
The probabilities, Pa{X^ Y), can be obtained in a similar way 
as in B-ICSK systems, and only Pa{0, 0) is shown here as an 
example: 

Pa(0, 0) = P5(0, 0, 0) + P5(l, 0, 0) + P5(2, 0, 0) + P5(3, 0, 0) 
= [\) [P(A^n<ri)+3P(7Vp + 7Vn<ri)" 



3-3Q 



ri - n{p2 -pi) 



y^n[p2{l -P2) -Pi)] +cr' 

For other expressions, please see Appendix |A| 



C. Molecular-Type-Based 

When different types of molecules represent different sym- 
bols, the technique is known as molecular- type-based modula- 
tion, referred to as molecule shift keying (MoSK) [ [T7| . Unlike 
the work in 1 17], we use here a set of isomers. We thus name 
it isomer-based MoSK (IMoSK). The IMoSK system requires 
only one threshold to be detected by a receiver nanomachine. 



which makes it simpler than the CSK system. For systematic 
analysis, we can choose one among the several aforementioned 
isomer sets; a modulation order can be determined by the sets 
used. For example, hexoses have a modulation order of up 
to 32, trioses of 4. For simplicity, in the first we apply, as 
was done in 1 17J , the AWGN model. The mutarotation effect 
explained in Section III-A is also considered in the second 
part. 

1) B-IMoSK-AWGN: Only AWGN considered. 



Pa{a, a) AWGN = Pb{(^, Q^, c^) + PbiP, c^, c^) 

= (2) [P{Np^N,^Nr,>r)P{Nr^<r) 
^P{N,^Nn>r)P{Np^Nn<r)\, 

Pa{(^, P)aWGN = Q^, P) + Ph{P, Q^, P) 

= (2) [P{Nn>T)P{N^^N,^Nn<T) 

+ P(7Vp + TVn > r)P(7Vc + A^n < r)] , 

Pa{P, ol)awgn = Ph{(^, (y) + Ph{P, P, (y) 

= [I) ' [P{Np + TVn > r)P(7Vc + A^n < r) 
+ P{Nr, > r)P{Np + AT, + AT, < r)l , 



)awgn Ph{0L,P, 



.p. 



P{N, 



P{N, + AT, > r)P(Ar^ + AT, < r) 
f AT, + A^, > r)P(Ar, < r) . 



2) B-IMoSK-Mutarotation Considered (a-D-glucopyranose 
<^ P-D-glucopyranose): When a- and P- I)-glucopyranose 
are chosen for the system, possibilities arise, due to mu- 
tarotation effect, of incorrect decoding, such as a- sent, P- 
received, or p- sent, a- received. Thus, we propose using the 
following expressions to calculate error probabilities, Pa{(y^ P) 
and Pa{P-, <^), considering the mutarotation process. The a- or 
P- form sent varies its number with time | [2T| , and the number 
can be calculated by observing the specific optical rotation. 
Rta or Rtjs (observed specific optical rotation at time t) minus 
Req (at equilibrium) divided by R^ ov Rfs (at time 0) minus 
Req has a linear relationship with time (Tg) at about 36.5° C, 
body temperature. Rta and Pt/3 values are calculated by R^ 
and Rj3, and Req can be found in pT| . The number of a- and 
P forms existing after time T5 is obtained as shown below. 
If the number of P- form exceeds the threshold value after 
Tg when a- sent, - value is added to the error term. Thus, 
Pa (a, P) can be calculated as 



Rh 



R. 



eq 



Pa Pe 



0.99 

3600^ 



Rf, 



Tig ^ n 
{Rtg - R^)n 
Pa — Pb 



ni3 



n — Tin 



If n(3 > r, Pa{a,P) = Pa{(y,P)AWGN 



n/3 
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TABLE I 

Comparisons of Modulation Techniques. 





Rate 


Tx Complexity 


Rx Complexity 


Distortion Sensitivity 


IMoSK 


good 


multiple kinds of messenger molecules 


multiple receptors 


robust 


ICSK 


moderate 


one kind of messenger molecules 


one receptor 


highly sensitive 


IRSK 


moderate 


two kinds of messenger molecules 


two receptors 


less sensitive 



Also, Pa{f3^a) can be calculated as 



Re 



R 



t(3 



0.99 



Rp 



:t. 



1, 



R 



t(3 



-R(3 3600^ 
_ UaRa + ni3Ri3 _ {u - n^)Ra + 

{Rt(3 - Ra)n 



R(3 — Ra 

If Ua > r, Paif^.a) 



n — rif^. 



Pa{P, Oi)aWGN 



where, Req = 52.7°, Ra=n22°, and R^=1^J° jlTJ. The 
parameters and n/3 indicate the number of a- form and 
/3- form molecules, respectively, and n is the number of total 
molecules transmitted. Optical specific rotation of the chemical 
compound is defined as the observed angle of optical rotation 
when plane-polarized light passes through the solution (i.e., 
I)-glucopyranoses). |^ 

3) 32-IMoSK: When using hexoses, the system has a 
maximum modulation order of 32. Probabilities of X sent and 
Y received for all X and Y values are obtained similarly. 

If X=Y (i.e., if X sent, X received.), 

32 

Pa{X,Y) = Y,Pb{Z,X,Y) 



PiNp + Ne + N„> t)P{N^ < t) 



31 (P(iVe + iV„ > T)P{Np + Nn< r)) 



If X^Y (i.e., if X sent, Y received), 

32 

Pa{X,Y) = Y.P'>{Z,X,Y) 



P{Np + iV„ > r)P(iVe + Nn<T) 



31 (P(iv„ > T)P{Np + + iV„ < r)) 



D. Molecular-Ratio-Based 

In addition to the concentration-based and the molecular- 
type-based modulations, this paper also suggests a new modu- 
lation technique as ratio-based modulation that can be referred 
to as isomer-based ratio shift keying (IRSK). IRSK encodes 
the information based on the ratios of messenger molecules. 

^It could be measured by a polarimeter, and there is a linear relationship 
between the observed rotation and the concentration of the compound (22|. 



Deploying this technique yields several benefits. First, it can 
have a high, theoretically infinite, modulation order. Moreover, 
we need only two types of molecules in the simplest system 
though it is possible to deploy many types of molecules for 
more complex systems. Even though it seems quite similar to 
the ICSK technique, the ratio-based technique can be simpler 
at the detection site since only ratios, not the concentration, 
need be detected. It is also more robust than the ICSK 
system when a channel is distorted (it will be specified in 
Section |III-E| ). As can be seen in Section |II-A[ the number 
of messenger molecules from the current and the previous 
symbols can be represented as a binomial distribution, and 
eventually, a normal distribution. Since we use two types of 
molecules in one symbol duration for a simple IRSK system, 
the notations can be modified as follows: 

Nci A/'(npii,npii(l -pn)), 

Npi ^ J^{npi2,npi2{l -P12)) - A/'(npii,npii(l -pn)), 
Nc2 ^ ^{np2i,np2i{l -P21)), 

Np2 J\fnp22,np22{'^ -P22)) - A/'(np2i,np2i(l - P21)) 

where, Nd and A^^i are the number of the first molecule from 
the current and the previous symbols received in the current 
symbol duration, respectively, and Nd and Npi are of the 
second molecule. Also, pij represents the hitting probability 
of zth molecule for jTg. 

There is one tip for picking up the proper kinds of messen- 
ger molecules to deploy the IRSK system. As mentioned in 
Section |III-C2[ mutarotation can, during propagation, occur in 
pairs of a and {3 forms. Thus, even though the effect is not 
significant, it is better not to choose the pair as the messenger 
molecules. A binary IRSK (B-IRSK) system is the same as a 
B-MoSK in that two symbols can be represented as the ratio 
1:0, or 0:1 of two types of molecules. Thus, we here analyze 
a Quadrature-IRSK (Q-IRSK) system, which is addressed in 
Appendix |B] and compare the result with those of ICSK and 
IMoSK systems. 

E. Comparisons 



The three techniques in Sections III-B III-C and III-D 



have pros and cons, and one must be chosen depending on 
several conditions. The comparisons are summarized here. 
First, achievable rates can be better for IMoSK systems unless 
it is binary case. This is because ICSK and IRSK systems 
have more thresholds, which means shorter minimum distance 
(i.e., concentration difference) between symbols. The IMoSK 
system, however, has to generate and detect multiple kinds of 
messenger molecules whereas the ICSK and IRSK systems 
need generate only one or two. Thus, transmitter and receiver 
complexities are higher in IMoSK systems. Lastly, channel 
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TABLE II 
Simulation Parameters. 



Parameters 


Values 


Ts |16 | 


5.9sec 


d 1 17 1 


16 jim 


Phit for Ts 


0.6097 


Phit for 2Ts 


0.7208 


Radius of the hexoses | 23 1 


0.38 nm 


D of hexoses jsj 


597.25 iim^/sec 


^-f^hexose 


nil kJ/mol 


Viscosity of the water 


O.OOlfc^f/sec • m 


Temperature (body temperature) 


36.5 °C = 310K 



i) Previous symbol = 'o' 



2) Previous symbol = 





Fig. 8. Binary symmetric channels of the binary channel model depending on 
the previous symbol (Z). The X and Y denote the transmitted and the received 
symbol, respectively, and Pb(Z, X, Y) is the probability of Z previously sent, 
X currently sent, and Y received. 



condition should also be considered to apply the proper 
modulation technique. That is, when the channel is distorted 
by any factor, the three systems are affected differently. In the 
worst case, there might exist some materials in the medium 
(i.e., channel) influencing the number of molecules named a. 
Then, the IMoSK system deploying molecules a, b, c, and d 
has an additional error probability of |, ICSK deploying only 
molecule a has an additional error probability of 1, and IRSK 

deploying the molecules a and b has the error probability of 

1 

2 • 

As specified above, each system has its advantages and 
disadvantages, and the most appropriate technique can be 
found given different situations. It will be numerically proven 
in Section IV Please see Table |I] for the summary. 



IV. Numerical Results 

Assume that all the hexoses have the same physical proper- 
ties: size, diffusion coefficient, and the enthalpy of formation. 
Phit value is calculated by the same numerical calculation used 
in [161 , it is approximated for the hexoses by assuming 
it to be proportional to the diffusion coefficient. Therefore, 
we have the hitting probabilities shown in Table |ll] and apply 
these to the proposed techniques explained in Section |Ill| Here, 
we define the achievable rate R that maximizes the mutual 
information I{X;Y) as follows: 



R 



X Y 

max/(X;Y') 



(7) 



E 



at 



0.2 



! 

: 




1 1 1 1 






i 1 1 ll>H-i 1 1 1 1 h 


- ; 




: 


jf^ gain 7 : 










— 1 — B ICSK (Proposed) 
^^^B-CSK (Conventional) 



10 



15 20 
SNR (dB) 



25 



30 



35 



Fig. 9. Achievable rate comparisons of the conventional insulin based CSK 
system and the proposed ICSK system using one kind of aldohexose isomers. 



E 

S 0.8 
s. 

^ 0.6 
■= 0.4 

I 

■g 0.21- 



I 


1 1 1 1 












^gain V 






i 




— 1 — B IMoSK (Proposed) 

B-MoSK (Conventional) 



10 



15 20 
SNR (dB) 



25 



30 



35 



Fig. 10. Achievable rate comparisons of the conventional insulin based 
MoSK system and the proposed IMoSK system using two kinds of aldohexose 
isomers. 



where, P{X) and P{Y) are the probabilities of events X (the 
transmitted symbol) and Y (the received symbol). 

From (|7]), we obtain the achievable rates under various 
signal-to-noise ratio (SNR) values. The SNR, we use here, 
is defined as the ratio of the received signal power/energy to 
the noise power/energy. The signal power is calculated from 
the number of messenger molecules transmitted, the hitting 
probability, and the energy model derived in Section |Il| Thus, 
the received number of messenger molecules is calculated 
using the transmitted number of molecules and the hitting 
probability, and the number is converted to power/energy 
unit by energy model. For binary channel analysis, binary 
symmetric channel (BSC) is used as shown in Fig. [5] 

Figs. |9] and [To] compare the achievable rates of the proposed 
method using hexoses (i.e., ICSK and IMoSK) with the 
conventional insulin-based method (i.e., CSK and MoSK) in 
binary systems. It is remarkable that, for each system, we 
obtain about 7 dB and 7.5 dB of SNR gain. This is due mainly 
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Fig. 11. Achievable rate comparisons of 32-IMoSK and B-IMoSK. The Fig. 13. Achievable rate of binary systems using hexoses as messenger 

32 kinds of isomers of aldohexose used for 32-IMoSK, and two kinds for molecules. 

B-IMoSK. 




to the size of the proposed messenger molecules being much 
smaller than that of insulin. Hence, the transmit energy of the 
proposed method is much less than that of the conventional 
method. Also, the diffusion coefficient is much larger for 
hexoses, which means faster propagation. 

In addition, the B-IMoSK system shows a better data rate 
performance than B-ICSK. Fig. [TT] shows the achievable rates 
of 32-IMoSK and B-IMoSK using hexoses as messenger 
molecules. As can be seen from Fig.[TT] the 32-IMoSK system 
has a maximum data rate of 5 (bits per symbol), and that of 
the B-IMoSK is 1 (bit per symbol). In Fig.[T2j we compare the 
achievable rates using trioses and hexoses. Obviously, trioses 
achieve higher SNR gain than hexoses due to their smaller 
sizes. Trioses, however, have the data transmission limit of 
one bit per symbol. From this result, we can conclude that 
trioses can be selected for a low data rate system with a higher 
transmission reliability and hexoses for a high data rate system. 



Figs. 13 and 14 compare the achievable rates of ICSK, 



IMoSK, and IRSK systems. In the binary case, the IRSK 
system is conceptually the same as the IMoSK while the 
rate goes the same way with the ICSK in the quadrature 
system. Fig. [15] especially shows the comparison in terms of 
modulation orders given SNR of 10 dB. The rate is not only 
higher in the IMoSK system, but the gap also increases slightly 
with modulation order. This can be caused by decreasing 
values between thresholds as modulation order increases. 

V. Conclusions 

To make nano communication feasible in practice, this 
work proposed novel modulation techniques using isomers as 
messenger molecules. We first introduced energy and channel 
models for our system. Next, several modulation methods were 
proposed so as to support up to five bits per symbol. We 
also compared the achievable rate performance with existing 
modulation methods (concentration-based and molecular-type- 
based), and suggested another modulation technique (ratio- 
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4 5 6 
Modulation Order 



Fig. 15. Achievable rate of IMoSK and ICSK systems in terms of modulation 
order when SNR is 10 dB. 



based). This work differs from prior work in that it proposes 
practical messenger molecules and provides guidelines for 
selecting from among several possible candidates. For the 
future work, we will consider deploying multiple sets of 
messenger molecules and collisions among them. Also, a 
velocity model as well as diffusion can be added to further 
increase the data transmission rate. 
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Appendix 

A. Q-ICSK 

Some probabilities, Pa{X^Y), for Q-ICSK systems are 
given as follows: 

Pa(0, 1) = P5(0, 0, 1) + P5(l, 0, 1) + P5(2, 0, 1) + P5(3, 0, 1) 
= [^) [P{ri <Nn< T2) + 3P(ri < ^ < T2) 

- Q (-) - 3Q I 

^'^^ Wn[p2(l -P2) +pi(l -pi)] +0-2 



Pa{0, 3) = Pfc(0, 0, 3) + Pb(l, 0, 3) + Pfc(2, 0, 3) + Pb(3, 0, 3) 

= (-] [P(iV„ > T3) + P{Np + Ar„ > T3) 



A/n[p2(l -P2) -Pi)] 

Pa(l, 0) = Pb{0, 1, 0) + Pfc(l, 1, 0) + Pb{2, 1, 0) + Pfc(3, 1, 0) 
= (-J \P{Ne + Nr,< n) + 5P{Np + N,+N„< n) 

1-Q 



3-3Q 



Ti — npi 
i/npi(l -pi) + 0-2 

Ti - np2 



V«[P2(l-P2)+2pi(l-pi)]+a2 



Pa(l, 1) = Pb{0, 1, 1) + P6(l, 1, 1) + Pb{2, 1, 1) + P6(3, 1, 1) 
= [^) [P(ti < A^c + A^n < T2) + 3P(ti <Np + N, + Nn< T2) 

+ 3| Q 



r2 - npi 



i/npi(l -pi) + 0-2^ 
Ti - np2 



i/npi(l -pi) + 0-2 



yn[p2(l-P2)+2pi(l-pi)]+(T2^ 

7-2 - np2 \ ] 

v/nb2(l-P2) + 2pi(l-pi)]+a2 ] j 



Pa(l, 2) = Pb(0, 1, 2) + Pb(l, 1, 2) + Pb(2, 1, 2) + Pb(3, 1, 2) 
= j [P(t2 < A^c + A^n < T3) + 3P(t2 < iVp + at, + Ar„ < Ts) 



3| Q 



T2 - np2 



y^npi{l -pi) +cr2 



T2 - n(p2 -Pl) 



Vn[p2(l - P2) + 2pi (1 - pi)] + a2 ^ 
Vn[p2(l -P2) + 2pi(l - Pl)] + (72 ] j 



P,(l, 3) = P5(0, 1, 3) + P5(l, 1, 3) + P5(2, 1, 3) + P5(3, 1, 3) 

' /l\'r 1 
= (^) [P{Nc + A^n > rs) + 3P{Np N, ^ > rs) 

-0 



T3 - npi 



P„(0, 2) = Pfc(0, 0, 2) + Pb{l, 0, 2) + Pfc(2, 0, 2) + Pfc(3, 0, 2) 
= j [P(t2 < A^„ < r3) + 3P(r2 < Np + < T3) 
2 ' T2-n{p2-pi) 



3Q 



A/npi(l -pi) +(t2^ 
r3 - np2 



Vn[p2(l-P2) + 2pi(l-pi)]+. 



g(^) -3Q 



\/^[p2(l -P2) -Pl)] +cr^ 

rs - n{p2 -pi) 



\/n[p2{l -P2) -pi)] ^(y^ 



jB. Q-IRSK 

The probabilities, Paisx^sy), for Q-IRSK systems are 
given, and here, 5i,52,ss, and 54 denote four different sym- 
bols (i.e., symbol set, S = {si, 52, ss, 54}). 
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[20] P. M. Dewick, Essentials of Organic Chemistry: For Students of Phar- 

-j 2 macy, Medicinal Chemistry, and Biological Chemistry. Wiley, 2006. 

P (sv Sv) = I — I \ ^ Phisy Sv Sv) ^^^^ ^' Carruthers and D. L. Melchior, "Transport of a- and /3-D-glucose by 

a\ X, YJ y^J b\ X, YJ .^^^^^ human red cell," Biochemistry, vol. 24, no. 15, pp. 4244-4250, 

Ph{sz, Sx.Sy) = P(rY_i < aiNpi + jSiNd + N^i < ry) [22] J. R. Mohrig, C. N. Hammond, and P. F. Schatz, Techniques in Organic 

Chemistry. Freeman and Company, 2010. 
X P(r_Y+4 < a2Np2 + f^2Nc2 + Nn2 < T_y+5^23] S. G. Schultz and A. K. Solomon, "Determination of the effective 

hydrodynamic radii of small molecules by viscometry," /. Gen. Physiol. , 

where, tq and indicate -oo and oo, respectively. vol. 44, no. 6, pp. 1189-1199, July 1961. 

If sz = si, then ai = 0, and a2 =1 . 

If Sz = S2 or 53, then ai = 1, and a2 = 1. 
If Sz = S4, then ai = I, and 0^2 = 0. 
If sx = si, then /Si = 0, and (32 =1 . 
If sx = S2 or 53, then /3i = 1, and (32 = 1. 
If Sx = 54, then /3i = 1, and f32 = 0. 
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